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Abstract: We present a new approach to combine multiple NLO parton-level calculations 
matched to parton showers into a single inclusive event sample. The method 
provides a description of hard multi-jet configurations at next-to leading order in 
the perturbative expansion of QCD, and it is supplemented with the all-orders 
resummed modelling of jet fragmentation provided by the parton shower. The 
formal accuracy of this technique is discussed in detail, invoking the example of 
electron-positron annihilation into hadrons. We focus on the effect of renormalisa- 
tion scale variations in particular. Comparison with experimental data from LEP 
underlines that this novel formalism describes data with a theoretical accuracy 
that has hitherto not been achieved in standard Monte Carlo event generators. 



1 Introduction 

During the past decade, Monte-Carlo methods for simulating hadronic final states in collider experiments 
have improved continuously. Multi-purpose event generators incorporating the most recent higher-order 
perturbative QCD calculations have thus emerged, making them available to phenomenology and experiment 
alike. This has far-reaching consequences for both precision physics and searches for new phenomena. Key 
to the developments has been the steady progress in understanding the interplay of real and virtual higher- 
order QCD corrections on one hand and of rcsummation techniques like parton-shower algorithms on the 
other hand. The construction and development of simulation tools for QCD processes has become one of 
the central activities of research in collider phenomenology. 

This publication discusses an extension to the established techniques of multi-jet merging and next-to- 
leading order matrix-element matching. Existing multi-jet merging methods (MePs) combine leading-order 
matrix elements of varying final-state multiplicity with parton showers. They were pioneered in [l}|4] and 
further matured in pUS]. The key advantage of these methods is the possibility to describe arbitrarily 



complex final states at leading order in the strong coupling, while providing fully inclusive event samples 
with resummation effects taken into account. They have therefore become standard analysis tools for collider 
experiments. However, they lack the precision of a full next-to-leading order perturbative calculation. This 
is remedied by next-to- leading order matrix-element matching methods (Mc@Nlo), which combine NLO 
QCD calculations of fixed jet multiplicity with parton showers. They were introduced in [9j|To] and have 
recently been automated in various programs |11||12| . Their main advantage lies in the excellent description 
of well-defined, inclusive final states. Using the MENLOPS technique (13 14 , it is possible to make these 
results exclusive and combine them with higher-multiplicity leading-order predictions in order to recover the 
virtues of MePs methods. 

The aim of this paper is to further improve upon the existing algorithms and construct a consistent, process- 
independent merging method for matched NLO predictions with varying jet multiplicity. Pictorially speak- 
ing, we intend to replace the leading-order matrix elements of the original MePs approach with corresponding 
ones at next-to-leading order. This is achieved by combining various Mc@NLO event samples and accounting 
for potential double counting by means of a modified truncated parton shower [5j[To] . Ultimately, we intend 
to maintain the fixed-order accuracy of the matrix elements, but also to preserve the logarithmic accuracy of 
the parton shower. The new method discussed here goes well beyond the scope of the MENLOPS technique. 

In the framework of this paper the formalism is specified for a multi-jet merging at NLO accuracy for 
-annihilations into hadrons, building on the existing implementations of MePs [5] and MC@NLO 
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techniques in the SHERPA event generator 15 16 . In the present paper, however, we will assume that the 
evolution parameter of the parton shower is defined in such a way, equivalent to the measure of hardness of a 
parton splitting, that effects due to a mismatch of these two quantities can be neglected. In other words we 
will neglect effects that arise from allowed emissions generated by truncated parton showers. An algorithm 
with the same goals and a similar setup for the parton shower has been detailed, also for e + e _ -annihilations 



into hadrons, in 17 . We note in passing that a more recent attempt to address this problem for hadronic 



initial states 18 could not be proven to maintain the NLO accuracy. It will therefore not be discussed. 



Instead, the reader is referred to a parallel publication, 19 , where the new formalism introduced here is 
discussed for the case of hadronic initial states, including all effects of truncated parton showers. As example 
there, highlighting the improvements due to this method over previous results, the production of W bosons 
in conjunction with jets, is considered. 

The outline of the present paper is as follows: Section [2] discusses the MePs algorithm for matrix-element 
merging at leading order, and the Mc@NLO method for NLO matching as implemented in SHERPA. As an 
intermediate step, the implementation of the MENLOPS idea for Mc@NLO core processes is presented. With 
the notation thus established, the new merging method at next-to leading order, MePs@NL0, is introduced 
in Sec.[3| The renormalisation scale dependence of the result is discussed in some detail. Sec. [4] is devoted to 
details of the Monte-Carlo implementation with particular emphasis on colour-correctness. Example results 
for the case of electron-positron annihilation into hadrons are shown in Sec. [5j including the impact of scale 
variations and of varying the number of jets described by NLO matrix element calculations. Sec. [6] presents 
our conclusions. 



2 Brief review of merging and matching techniques 

In this section, existing matrix-element parton-shower merging and matching methods are briefly reviewed, 



using the notation of 12 . As already stated in the introduction, the effects of allowed emissions generated 
by truncated showers are ignored, which improves the readability of this publication, allowing to focus on 
the structure of the result. For a full algorithmic solution, we refer to the parallel publication, in |19| . Our 
approach is justified by the choice of transverse momentum as evolution variable in the parton shower used 
for this publication. 

In the context of merging, we define a jet criterion Q n , which typically denotes the minimal value of some 
relative transverse momentum present in the n- parton phase-space configuration <!>„. Correspondingly, Q cut 
denotes a jet-defining cut value, called the merging scale, such that for n-jet events the condition Q n > Q cu t 
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is fulfilled 

Formally, the quantity of interest is the expectation value (O) of an arbitrary, infrared-safe observable O, 
evaluated by taking the average over sufficiently many points in an n-particle phase-space, 

The methods reviewed here, as well as our newly proposed technique, have the following aims 

• Multi-jet merging techniques 

For configurations with n jets, the respective fixed-order accuracy of (O) inherent to the parton-level 
result should be maintained. More precisely, for leading-order merging (MePs), jet observables for n 
jets above the merging scale Q C ut should be determined at leading-order accuracy. For next-to-leading 
order merging (MePs@Nlo) they should be given at NLO accuracy. For configurations below Q C ut, 
the MePs accuracy will be that of the parton shower, while for MePs@NL0 leading-order accuracy 
is envisaged. At the same time we require that the logarithmic accuracy of the shower be either 
maintained or improved in the region above Q cu t. 

• NLO matching methods 

For processes leading to n-parton final states at leading order all n-particle inclusive observables, and in 
particular the total cross section, are expected to reproduce the fixed order NLO results. At the same 
time, all n+l-particle observables are expected to be given at leading order accuracy, while higher-order 
emissions should still be described by the leading logarithmic approximation of the parton shower. 



2.1 Leading-order merging - MePs 

In the context of the leading-order merging method proposed in [5], the following quantities are introduced: 

• Squared leading-order (Born) matrix elements, B„($„), for n outgoing particles, summed (averaged) 
over final state (initial state) spins and colours and including symmetry and flux factors. 

• Sudakov form factors of the parton shower, given by 

(2-2) 



A r ( l PS )(<,i')=cxp|-y'd$ 1 K„( 



K„ denotes the sum of all splitting kernels for the n-body final state. The one-particle phase-space 
element for a splitting, d^i, is parametrised as 

d$i = dtdzd(/) J(t,z,<t>) , (2.3) 

where t is the ordering variable, z is the splitting variable and <fi is the azimuthal angle. J(t, z, </>) is 
the appropriate Jacobian factor. The ordering variable is usually taken to fulfil t oc as t — > 0. 

The resummation scale «q, which defines an upper limit of parton evolution in terms of the shower 
evolution variable. t c is an infrared regulator of the order of Aqcd marking the transition into the 
non-perturbative region. 



The expectation value of an arbitrary, infrared-finite observable O to 0(a s ) has been computed in 14 . It 
is derived from the following expression: 

2 

r 

(O) = y"d$„B„($„) A( l PS )(i c ,^)0($„) + y"d$ 1 K n ($ 1 )A( PS )(i„ +1 ,^)e(Q cut -Q n+1 )0($„ +1 ) 
+ J d$„ +1 B n+1 ($„ +1 )A( PS )(t n+1 ,^)0(Q„ +1 -Q cut ) 0($„+i) , 



1 The jet criterion Q applied here has been slightly modified compared to [H], in order to reflect the fact that no unique 
parton flavour can be assigned at the next-to-leading order. For any pair of final-state partons i and j we define 

Qij = ZPiPj min 2 where C* = 7 — p ' Pk . (2.1) 

The spectator index k runs over all possible coloured particles. 
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(2.4) 



where (9($ m ) is the observable evaluated for an m-parton final state. The term on the second line yields 
leading-order accuracy for any n + 1-particle observable in the region Q n +\ > Q CU f Leading-order accuracy 



for observables sensitive to $„ is guaranteed by the fact that Eq. (2.4) can be rewritten as 



(O) 



d$„ B n ($„) 



A( l PS )(t c , M |)0($„) + J d*i ^(iilAf'^.^JO^) 



(2.5) 



d$ 



n+l 



B n +i(*„+i) - B n ($„) K„($„+i) 1 A^ PS) Mq) 9(Q«+i - Qcut) 0($„ +1 ) , 



where the first line is the 0(a e ) parton-shower result [20] and independent of Q cu t. The additional terms 
on the second line incorporate possible sub-leading colour single logarithms as well as power corrections. 
The size of these corrections determines the potential discontinuity in (O) at Q cu t ■ It can be large if Q cu t is 
cither far from the collinear limit or sub- leading colour single logarithms are important. Sub-leading colour 
configurations, however, can be included in a systematic manner, as was detailed in [12] . 



An important feature of Eq. (2.4 1 is that it can be iterated to incorporate higher-multiplicity leading-order 
matrix elements into the prediction. By replacing n — > n + 1, all properties of the algorithm remain the 
same. In order to obtain this property when dealing with next-to-leading order matrix elements, a slight 
modification is necessary, which will be described in Sec. [3] 

2.2 Next-to-leading order matching - McONLO 

In the Mc@NLO matching method the following additional quantities are needed: 

• Squared real-emission matrix elements, R„($„+i), for n-particle processes, summed (averaged) over 
final state (initial state) spins and colours and including symmetry and flux factors. Note that 
Rn($n+l) = B„+i($ n+ i). 

>(A) 



The NLO-weighted Born differential cross section B„ , defined as 
- B n ($ n ) + V n ($ n ) + l! s) ($«) 



d$i 



Di A >($„ +1 ) e(/4 - t n+1 ) - Di s )($„ +1 ) 



(2.6) 



(S) 

Here, V„ i s th e Born-contracted one-loop amplitude, and In is the sum of integrated subtraction 
terms, cf. 
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while D { n ] are the corresponding real subtraction terms. In contrast, D,^ are the 
Mc@NLO evolution kernels multiplied by Born matrix elements. Both functions can be decomposed 
in terms of dipole contributions, D = J^ij k ^ij,k, where each dipole encodes exactly one singular 
Further, each dipole has a corresponding phase space factorisation d$„ +1 = d$„ d$ 
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region 

and t n +i = t($„-|_i) is defined in terms of Eq. (2.3 1 in each of these dipole phase spaces 



ij,k 



• The hard remainder function 

H! A >($„ +1 ) = R n (<f n+1 ) - D,[ A )($„ +1 ) 9(/x| - t n+1 ) . 
with t ]v+i = t(<f> n+ i) defined as above. 

• The McQNLO Sudakov form factor 

d4 a) ($„,$i)1 



(2.7) 



a! a) (m') 



exp 



d$i 



B(*„) 



(2- 



(A) (PS) 

Note that A„ implicitly depends on <!>„, while the original Sudakov form factor A„ does not. This 
is a consequence of the fact that the two Sudakov form factors differ by their treatment of colour and 



spin correlations and it is needed to maintain full NLO accuracy. For a detailed discussion, see 12 
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The expectation value of an arbitrary infrared safe observable O to 0(a s ) is then given by [9] 

r ^ Q ( A ) 

(O) = J d$„B( A )($„) A< A >(i c ,^)0(<I> n ) + J d$! D ' B ^^) $l) Al A) (Wi^Q) 0($n+l) 
+ y d$„ +1 H^($ n+1 ) 0($„ +1 ) . 



(2.9) 



Events generated according to the first line are referred to as standard, or S-events, while events generated 
according to the second line, the hard remainder, correspondingly are dubbed H-events [9 12 . Note that 
the square bracket in the first line integrates to one, reflecting the probabilistic nature of the Sudakov 



form factor. This, together with equations (2.6 1 and (2.7), implies that the total cross section reproduces the 



exact NLO result. Correspondingly, an Mc@NLO prediction is next-to-leading order accurate for observables 
sensitive to the Born phase-space variables $„, and leading-order accurate for observables sensitive to $ n +i- 
It should be stressed at this point that this statement is true for arbitrary processes with arbitrary colour 
structures only if the first emission in the parton shower is generated by terms D„ ' . This modification of the 
parton shower, as achieved in the SHERPA implementation, allows to guarantee fully colour accurate results. 
Such terms have neither been included in the original implementation [9] nor in the more recent version 
presented in 11 . It also is worth mentioning that, in contrast to the MePs method, leading-order accuracy 



is maintained throughout the n + 1-particle phase space, but it cannot be extended to higher parton or jet 
multiplicity. This extension will be the topic of Sec. |2.3| 



The second line of the equation above, (2.9), may be supplemented with a Sudakov form factor, 



(0)= |d$„B( A )($„) Ai A \t c ,4)0($ n ) + J dS x \ ( ^ 1) a 1 A) (Wi,^)0('I' 
+ J d<E>„ +1 Hi A >($„ +1 ) Ai PS )(t n+1 ,/4) 0(* n+1 ) . 



(2.10) 



While this does not alter the fixed order accuracy of the method, it improves its logarithmic behaviour, 
which can now be seen to be identical to that of the underlying parton showeij^] This alternative form has 
also been implemented in SHERPA. 



2.3 Combining NLO matching and LO merging - MENLOPS 



NLO-matched predictions as described in Sec. 2.2 can easily be merged with higher-multiplicity event samples 
at leading order accuracy using the techniques described in Sec. |2.1| T he original algorithm, based on the 
POWHEG method 10 21 , was independently proposed in [l3| and [14] . In this publication we extend the 



method to Mc@NLO, which requires the introduction of the local K- factor 



Bl A) (^) 
B„(*„) 



1 - 



H„($ 



n+l, 



H n ($ M 



R„($„+i)y R„($ n+ i) 



(2.11) 



It is motivated by the behaviour of the underlying Mc@NLO event sample in terms of S- and H-events 9p2~ 



In both cases, however, the n-parton phase-space configuration in Eq. (2.11) is determined by backward 
clustering, as described in (HI. 

The expectation value of an arbitrary, infrared-finite observable to 0(a s ) in the MENLOPS method for 



2 Of course, if new Born-level configurations enter the hard-remainder, which cannot be accessed by the parton shower, they 
are not modified with such a Sudakov form factor. As an example for such a configuration consider the process uu — > W + sc 



at a hadron collider, which although part of the W+jets ensemble cannot be obtained from any W 
through parton evolution. 



1 parton configuration 
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Mc@NLO is given by 



J (2.12) 



(O) = J d$„B( A )($„) 

l ~ w (A) 

A,( A )(t c , / i|)0($„) + J d#l ""b^^^ A n A) (*n+l,Atq) e(Qcut - Qn+l) 

+ /" d$„ +1 Hi A )($„ +1 ) A( l PS )(i n+1 ,^)e(Q cut - Q„+i) 0($„+i) 

+ y"d$„ +1 fci A )($„ +1 )B„ +1 ($„ +1 )A( l PS )(t„ +1 ,^)8(Q n+1 -Q cut ) 0($ n+ i) . 

This prediction is next-to-leading order accurate for observables sensitive to $ n and leading-order accurate 
for observables sensitive to $ n +i- The key advantage compared to a pure NLO-matched prediction is 
that final states of higher jet multiplicity are treated as in the MePs approach. The improvement over 
results obtained by MePs methods is the next-to leading order accuracy of the inclusive cross section and 
of observables sensitive to $„. 

The method aims to maintain the full NLO-accuracy in the n-jet phase space and the LO-accuracy in the 
(n + l)-jet phase space, without upsetting the logarithmic accuracy of the parton shower. In order to see 
that this indeed is the case, equation (2.121 can be rephrased as follows: 



(O) = (O) 



MCONLO 



(Of 



(2.13) 



MCONLO 



given by (2.101, and thus showing the desired property. It thus remains to show that the 



with (O) 

correction term does not introduce unwanted terms of 0(a s ). Ignoring the obvious phase space arguments 
of the different matrix element contributions, it is given by 



(O) corr 



j d$„ +1 e(Q n+1 -Q cut )0(^ n+1 )A^ s \t n+uf j, 2 Q ) 



B 



(A) 



1 - 



H 



(A) 



tt(A) 
tin 



B„ 



+ i 



B n y B„ + i I B„ + i 
d^n+i Q(Q n +i -Qcut)0($ ti+1 )A( ps )(i n+1 ,^) 

x f 6 " A Wl A " A) ^ + ^q) M 



H (A) 



6(A) A (A) (t n ll 2 ) 

™ Ej( A ) \. L n+l>f*Qj 



B„ 



A n (t n _|_i, fjig) 



(2.14) 



Since D^ A ^ is of 0(a s ) and because the ratio of Sudakov form factor is at most of non- leading logarithmic 
order, 0(a s L), and non-leading in 1/N C , the overall contribution of this term is at most of 0(a 2 L), sub- 
stantiating the claim of accuracy. Higher jet multiplicities then exhibit the same accuracy as in the MePs 
approach. 



3 Merging at next-to leading order 

The previous section sets the scene to introduce a new prescription, which consistently merges multiple 
MC@NLO-matched event samples of increasing jet multiplicity. The method is constructed such that it is 
next-to-leading order accurate for observables that are sensitive to both <!>„ and $„+i 0(Q C ut — Q), while 
maintaining the logarithmic accuracy of Mc@NLO for observables sensitive to Q n +i- In other words, the 
goal is to describe every jet observable at next-to leading order in the strong coupling constant, including 
Sudakov suppression factors. 
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3.1 Definition of the MePsONlo technique 



We propose a method based on the following expression for the expectation value of an arbitrary infrared- 
finite observable O 



(O) 



D (A) 



J d$„ +1 H,[ A ) A,( l PS )(t n+1 , / i|) 6(Q cut - Q n+1 ) O n+1 

2 

- J d^ 1 K n jA^ s Ht n+1 , f i 2 Q )e(Q n+1 -Q cut )O n+1 



n+l \ D n+1 



n+l 
n+l 



(3.1) 



/D 
■t>n+l 

L to 

| d$ n+2 4 A \ a£_ { (* n+2) t n+1 ) A(f s )(t n+1 ,^ ) Q (Q n+1 - Q cut ) O. 



'n+2 + 



where again the obvious phase space arguments in the matrix element contributions and splitting kernels 
have been suppressed for better readability, and where they have been moved to subscripts in the observables. 
The dots indicate contributions from higher parton-level multiplicities, which are dealt with in an iterative 



procedure as detailed in Sec. 3.2 



It is easy to show that next-to-leading order accuracy is maintained for observables sensitive to & n +i a t 
Q > Qcut j where Q is the transverse momentum scale of the first emission, i.e. of parton n + l. Expanding 
the Sudakov form factor a! PS ^(£, /ig) in the third line to first order and combining it with the square bracket 
on the same line yields correction terms which are at most of 0(a 2 s ). 

In order to show the logarithmic accuracy of the procedure, Eq. (3.1 ) is rewritten as follows 

(O) = (0) MC@NL0 + (O) corr , (3.2) 



with ((^MCaNLO gj ven by (2.10|. Taking into account only n+l parton final states the correction term this 
time is given by[[] 



(oy 



y"d$„ +1 Q(Q n+1 - Qcut) 0„+i A[ p + s »(t cl Wi)A( ps '(t n+1) , 1 2 Q ) 



x < B 



(A) 
n+l 



B 



n+l 



(A) 



B 



^i+\(^c, t n+ i) ^ , A) ^ B„ Ai ( A) A n '(t n+ i,^Q) 



.(A)/ 



x <^ D 



J d$„ +1 6(Q„+i - Qcut) 0„+i &g B \t n+ i,v%) 



(3.3) 



1 - 



Bn A„ PS) (i n+ l,/4) 



— B n+ i 



1 - 



'r(A) 

B„+i 



A^ifte, Wi) 



Both terms in the curly brackets consist of one factor describing the emission of an extra particle, D n and 
B„ +1 . Those will eventually yield a contribution of 0(a s L 2 ). The factors multiplying these emission terms 
are at most of 0(a s L). However, these logarithms, if present, are due to sub-leading colour configurations 

3 Additional contributions are at most of 0(a^L 2 ) and thus do not impair the logarithmic or fixed order accuracy we intend 
to prove. 
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stemming from the difference between A^ A ) and A^ PS ^. The combination of virtual and real contributions in 

- (A) 

B,i does not induce any logarithms spoiling the accuracy of the parton shower. Thus the correction term 
does not impair the formal logarithmic accuracy of the parton shower. 

It is worth noting here that the algorithm detailed in 17 , while aiming at the same formal accuracy, follows a 
different construction paradigm. Rather than starting from the matrix elements, like the approach presented 
here, and matching the showers to them, its authors start from the parton shower and correct its emissions 
with higher order matrix elements. 



3.2 Iteration for multijet events 



Having shown, for the case of the first additional emission, how NLO- and the logarithmic accuracy of the 
shower are maintained, we now turn to the question how this can also be shown for the fcth additional jet. 
The first thing to be understood is that, in general, the observable O will have support in different sectors by 
different jet multiplicities. In the formalism outlined here this is reflected by the 0-functions involving the 
jet cut Q cu t and the scale Q of the softest emission of a given Borndike (n + fc)-jet configuration, in general 
given by Q n +k — Q($n+k)- For such a configuration, the respective expression for the (n + fc)-exclusive jet 
part of the observable, 



(0) C n X +k = E <°J 9 W"+fc - Qcut)e(Q cut - Qn+k+l)) 
j=n-{-k 



(3.4) 



is given by the suitably modified second part of Eq. (3.1 1, 
d$ n+k 6(Q n+fe -Q cut )B^_\ 

ti \ 

B„ + i 



n+k-l / 

n Af s »(t 1+1 ,i.) i 



(A) 



B 



n-\-k 



t i+ i 



D 



(A) 



Al^ k (t c ,t n +k) O n+ k + / d$i " +fc A^ fc (^„ + fc + i, t n+k )Q(Q cut — Qn+k+i) On+fc+i 

-Dn+fe 



to 



(3.5) 



/n+k 
d$ n +fc+l Q{Qn+k — Qcut) ©(Qcut — Qn+k+l) O n+k+ \Yi^ k ]^[ A^ ?S ^ (t i+ i , ti ) . 

i—n 

In order to see the formal accuracy of this expression, let us define an (n + fc)-jet inclusive expression of 
the observable, by dropping the second 0-function in (3.4 1. As before, it can be written as the sum of an 
Mc@NLOdike expression acting on the (n + fc)-parton Born configuration and a correction term, 

(oy:tk = (o)^r LO + wzii (3-6) 

where 

(A) 



(o)r + r LU = dt> n+k &(Q n+k -Q cut )B^ k 



' n+fc-l / 

n a^wo i 

i—n \ 



B 



n+k 



(A) 



B 



d$iK t 



'U+i 



A^ k {t c ,t n+k )O n+k + I d$i 



D 



(A) 

n+k a (A) 



B 



n+k 



A n+k (tn+k+l,t n +k) On+fc+i 



(3.7) 



« n+k 
/ d&n+k+l ©(Qn+fc — Qcut) O n+k +\ H n+fe 



The only difference with respect to the usual form of the Mc@NLO expression in (2.10) is the term in the 
second line which encodes a veto on emissions into the jet region from intermediate lines with its C(a s )-part 
subtracted. 
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At the relevant order in a s , this correction term reads 



(O) 



corr 
n-\-k 



/i t-prv-p J- 
d$ n+fc+1 Q(Q n+k+1 - Q cut ) O n+k+1 Y[ Af S) (t.- + i,t 



B 



n+fc+1 



B 



(A) 

n+fc+1 



B 



n+fc+1 



(A) n+fc-1 



1 - 




£ / d$!K 4 



A^ fe (i„+fe+i, t n+k ) 



^i+fc+lftc; ^n+fe+l) 
\l+2+l(*c 5 ^n+fc+l) 



(3.8) 



and the same reasoning already applied to Eq. (|3.3|) yields the desired result. For a more detailed discussion, 
including the effect of truncated showering, see |19| . 

The finding above shows that no terms appear due to the merging prescription that violate the logarithmic 
accuracy of the parton shower at and around Q cu t- To see this, it is sufficient to analyse the first emission 
off the (n + £;)-jct configuration over the full phase space. The second emission is, of course, completely 
determined by the parton shower and thus correct by definition. Also, clearly, the phase space for this first 
emission is confined to the region below Q cu t, therefore the behaviour above this scale is defined by the 
parton-level result with next higher multiplicity, the (n + k+ l)-jct configuration. By however extending the 
first emission above this cut and analysing the impact on O n+ k+i we show that the two regions match as 
smoothly as the logarithmic accuracy of the parton shower dictates. 



3.3 Renormalisation scale uncertainties 

The key aim of the MePs@NL0 approach presented here is to reduce the dependence of the merged prediction 
on the renormalisation scale ^t#, which is employed in the computation of the hard matrix elements. This 
scale has not been made explicit so far. 

Note that only the dependence on the renormalisation scale is reduced compared to the MePs method, while 
the dependence on the resummation scale, [Iq, remains the same. This is a direct consequence of the fact 
that the parton-shower evolution is not improved in our prescription, but only the accuracy of the hard 



matrix elements. The resummation scale dependence was analysed in great detail in 12 . 
Following the MePs strategy, the renormalisation scale should be determined by analogy of the leading-order 
matrix element with the respective parton shower branching history (5l. In next-to-leading order calculations, 
however, one needs a definition which is independent of the parton multiplicity. The same scale should be 
used in Born matrix elements and real-emission matrix elements if they have similar kinematics, and in 
particular when the additional parton of the real-emission correction becomes soft or collinear. This can be 



achieved if we define the renormalisation scale for a process of C(a") as 22 



a s (v 2 R y l = Y[a s (tf) , (3.9) 

i=i 

a procedure that has been used in LO merging for some time. Here, are the respective scales defined by 
analogy of the Born configuration with a parton-shower branching historjj^] 

The renormalisation scale uncertainty in the MePs@NL0 approach is then determined by varying /ir — > flu, 



while simultaneously correcting for the one-loop effects induced by a redefinition in Eq. (3.9). That is, the 
Born matrix element is multiplied by 

*s(&r (i - ^Ir 1 a E lQ g fj) ' ( 3 - 10 ) 

to generate the one-loop counter-term, while higher-order contributions remain the same. 



4 In the case of the real-emission correction and the corresponding dipole subtraction terms we consider the underlying Born 
configuration instead. The same scale definition is used in the parton shower and, consequently, in the Sudakov form factors. 
Of course, the nodal scales fii found in the backward clustering on the Born-like configuration of a single event then enter the 
truncated showering. 
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4 Monte-Carlo implementation 



In this section we describe the Monte Carlo implementation of the merging formula Eq. ( 3.1 1 in SttERPA. The 
techniques needed to combine leading-order matrix elements with parton showers are given elsewhere [5] . 

4.1 Generation of the parton-shower counterterm 

In addition, we now have to implement a method to generate the parton-shower counterterm on the third 



line of Eq. (3.1 ). Note that, by construction, this counterterm has the same functional form as the exponent 



of the Sudakov form factor A« (t, £*q)- We can therefore use the following algorithm: 

• Start from an n-parton configuration underlying the n + 1-parton event at scale /ig , 
and implement a truncated parton shower with lower cutoff scale t. 

• If no emission is produced, the original n + 1-parton configuration is retained. 

• If the first emission is generated at scale t' with Q > Q cu t , the event weight is multiplied 
by 1/k, where k — B^ l A h '' 1 ((f> n+ i)/B„ + i($„ + i). Evolution is restarted at t' . 

• All subsequent emissions are treated as in a standard truncated vetoed parton shower. 
Events will then be distributed as 

1 /""a 



Al PS >(^) 

= Ai PS >(i,/4) 



d$i 



K n ($r) 0(Q - Q cut ) A( ps )(t', tf Q ) A( ps )(M') 

2 

d$i K„($i) Q(Q - Q 

cut ) 



(4.1) 



This simple algorithm allows to identify the 0(a s ) counterterm with an omitted emission and to generate the 
correction term on-the-flight, much like the Sudakov form factor is computed in any parton-shower algorithm 
itself. 



4.2 Generation of the Mc@Nl_0 Sudakov form factor 

In this subsection we briefly recall an algorithm to compute Mc@NLO Sudakov form factors Ff2], as it is one 
of the basic ingredients to guarantee next-to-leading order accuracy of our MePs @NLO predictions. 
It is well known how to generate emissions according to Sudakov form factors with strictly negative exponent. 
In our implementation of Mc@NLO, however, we have to deal with potentially positive exponents, related to 
sublcading colour configurations. This leads to form factors larger than one, which cannot be interpreted in 



terms of no-branching probabilities and which are dealt with using a modified Sudakov veto algorithm 12 23 

Assume that f(t) is the sole splitting kernel in our parton shower, integrated over z and (p. The differential 
probability for generating a branching at scale t, when starting from an upper evolution scale t' is then given 
by 

T(t,t') = f(t) exp J - /' dtf(t) 1 . (4.2) 



The key point of the veto algorithm is, that even if the primitive of f(t) is unknown, one can still generate 
events according to V using an overestimate g(t) > f(t), if g(t) has a known integral. Firstly, a value t is 
generated as t = G _1 [G(t') + log#]. Secondly, the value is accepted with probability f(t)/g(t) 



24 



One can now introduce an additional overestimate h(t). The related weights are applied analytically rather 
than using a hit-or-miss method. They can thus be used to absorb the negative sign of the Mc@NLO kernels 
D« /B n . This leads to a correction factor for one accepted branching with m intermediate rejections of 

W (t, tl ,...,t m) = W) n mg{ti) _ m , (4.3) 
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where the ti run over intermediately rejected steps. Note that Eq. (4.3) can lead to negative weights, which 
reflect the fact that sub-leading colour configurations are taken into account and that the a-priori density 
h{t) might underestimate fit). 

In order to implement an evolution using the MC@NLO kernels D, ( l A) /B„ we need to identify the function 
/ above with the (2,</>)-integral of these kernels. A convenient choice of the function h will be the (z, <ft)- 
integral of the parton-shower evolution kernels K„. We are then free to choose the auxiliary function g on 
a point-by-point basis, but a convenient way is to define g = C f, where C is a constant larger than one. 
This guarantees that both acceptance and rejection terms are generated in sufficient abundance to reduce 
statistical fluctuations. 

The above method guarantees that all subleading colour single logarithmic corrections to B n are exponenti- 
ated. One can therefore guarantee the process-independent next-to-leading order accuracy of the Mc@NLO. 



5 Results 



15 16 



In this section results obtained with the MePs@NL0 method are presented for the case of e + e~ -annihilation 
into hadrons. The general-purpose event generator SHERPA sets the framework for this study 
Leading-order matrix elements are generated with AMEGIC++ 
subtraction 
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28 



25 and COMIX 26 . Automated dipole 



are employed to obtain parton-level events at 

The parton shower 



33 ; the related Mc@NLO generator has been 



and the Binoth-Les Houches interface 
next-to-leading order with virtual corrections provided by the BLACKHAT library |29f|32 
in SHERPA is based on Catani-Seymour dipole factorisation 
presented in [12] . In contrast to all other Mc@NLO implementations, no leading colour approximation is 
made and the full next-to-lcading order accuracy of the fixed-order result is therefore recovered throughout 
the phase space, cf. Sec. [4^2] The results presented here are at the hadron level. Note that the hadro- 
nisation model in SHERPA [34] has been tuned in conjunction with the parton shower and leading order 
matrix elements. It is therefore not surprising when deviations are found in observables that are sensitive 
to soft particle dynamics. In the future this will necessitate a new tune of the hadronisation based on the 
NLO-merging outlined here, rather than on the LO MEPS prescription that has been used so far in SHERPA. 



5.1 Choice of the merging scale 

Figure [T] shows the dependence of MePs@NL0 predictions for the Durham jet resolution on the merging 
scale Qcut- In order to match the customary notation we quote Y cut = {Q C ut/E cms ) 2 . All results were 
generated using 2-, 3- and 4-jet NLO parton-level calculations combined with 5- and 6-jet at leading order. 
The variation of results with Y cut in the region below and around Y cut is of the order of 10%, the predictions 
above the cut are remarkably stable and match the experimental data very well. Consequently, one should 
always choose the merging cut such that the analysis region is fully contained in the region covered by the 
NLO calculation of interest. 



5.2 Comparison of approaches and their perturbative uncertainties 

In this section we compare the renormalisation scale uncertainties between the MENLOPS and the MePs@NL0 
method. We choose jlp — C/j,r with C E {0.5, 1, 2} and set Y cut — 2. In the MePs@NL0 sample we generate 
2-, 3-, and 4-parton final states at NLO and 5- and 6-parton final states at LO. The MENLOPS sample only 
has the 2-parton final state at NLO and the remaining multiplicities up to 6 partons from tree-level matrix 
elements. Figures [2] to [8] show the respective scale variations as bands around the central prediction with 
C = 1. A significant reduction of the scale uncertainty is found for those observables, which are sensitive 
to the NLO parton-level results. This can be seen in particular in Fig. [2] where the 2 — » 3 and 3 — > 4-jet 
rates show significantly reduced uncertainties for larger y, while the 4—^5 and 5 — > 6-jet rates do not. 
Similar effects are observed in most event shape distributions in the hard region, for example in Fig. [3] for 
T — » 0.5. The reduction of the scale uncertainty in the moments of the event shape distributions in particular 
is more than impressive. It is also worth pointing out that the typical Sudakov shoulder at C = 0.75 in the 
C-parameter, which is notoriously difficult to describe in fixed-order calculations, now shows a remarkably 
smooth behaviour due to the successful interplay of the different multiplicity contributions. 

A final comment, concerning the evaluation of theory uncertainties by scale variations is in order here. 
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Durham jet resolution 3^2 (E CMS = 91.2 GeV) 



Durham jet resolution 4^3 (E CMS = 91.2 GeV) 




Figure 1: Experimental data from ALEPH 35 for the differential (n+1) — > n jet rates with n = {2, 3, 4, 5} 
(upper and lower panel, left to right) at the Z pole (i? c . m . — 91.2 GeV) are compared with 
MEPsONLO simulations with different values of the merging cut, F cut = lO^f 1 - 75 ' 2 ^ 2 - 25 !. To 
guide the eye, the merging cuts have been indicated with dotted lines in the same colour in the 
ratio plot. 
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Durham jet resolution 3^2 (E CMS = 91.2 GeV) 



Durham jet resolution 4 — > 3 (Ecms — 91.2 GeV) 




Figure 2: Perturbative uncertainties in MENLOPS and MePs@Nlo predictions of differential jet rates com- 



pared to data from ALEPH 35 . 



Clearly, there are two sources of perturbative uncertainties: the one analysed here, which stems from the 
matrix element. It is thus susceptible to variations of the renormalisation and, if present, the factorisation 
scale. In addition, changes in the value of a s , which we did not pursue here, or in parton distribution 
functions would have to be considered for a more complete assessment of such uncertainties. On the other 
hand, there are, of course, also uncertainties in the treatment of secondary emissions through the parton 
shower. There, in addition to the variations outlined above, one could also vary the parton shower starting 
scale, fiQ, which is equivalent to a variation of the corresponding resummation scale in analytical calculations. 
Obviously in regions that are dominated by the parton shower, such a variation would give a more sensible 
estimate of theory uncertainties than a variation of the scales in the matrix element, that we focused on here. 
As an example for this, consider the low-p± regime of the differential jet rates yij, — logj/jj — > oo. There the 
bands obtained from a scale variation in the matrix element regime are suspiciously small, and it is clear that 
a variation of the resummation scale would yield larger uncertainties. The same, obviously is true for the 
MePs@NL0 and the MENLOPS method, since in such regions both exhibit a comparable formal accuracy. A 
careful analysis of such effects, however, is beyond the focus of this paper, which discusses improvements of 
our ability to generate inclusive samples of events by increasing the formal accuracy of the matrix element 
part of the simulation. We therefore postpone this discussion to future work. 
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Thrust (Bcms = 91.2 GeV) 



Moments of 1 — T at 91 GeV 
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Figure 3: Pcrturbative uncertainties in MENLOPS and MePs@Nlo predictions of thrust. Compared are 
the measurements for the event shape from ALEPH [35] and its moments from OPAL 
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Figure 4: Pcrturbative uncertainties in MENLOPS and MePs@NL0 predictions of total jet/hemisphere 
broadening. Compared are the measurements from ALEPH [35] and OPAL 
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Figure 5: Pcrturbative uncertainties in MENLOPS and MePs@NL0 predictions of wide jet/hemisphere 
broadening. Compared are the measurements from ALEPH 35 and OPAL 36 . 
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C-Parameter (E CMS = 91.2 GeV) 



Moments of C at 91 GeV 




Figure 6: Perturbative uncertainties in MENLOPS and MePs@NL0 predictions of the C-parameter. Com- 
pared are the measurements from ALEPH [35] and OPAL 36 . 
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Figure 7: Perturbative uncertainties in MENLOPS and MePs@NL0 predictions of sphericity. Compared 



are the measurements from ALEPH 35 and OPAL 36 
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Figure 8: Four-jet angles using the Durham algorithm compared to data from OPAL 37 
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6 Conclusions 



In this paper we have introduced a method for a consistent multijet merging at NLO accuracy for the case 
of e + e~ -annihilations to hadrons. By explicit calculation, we have shown that our description maintains 
the higher order accuracy of the underlying matrix elements in their respective phase space range, while the 
logarithmic accuracy of the parton shower is respected. We have also analysed the impact of renormalisation 
scale variations in our new formalism. The results displayed here are exemplary for a far wider range of 
observables, which show a very good agreement between our simulation and data throughout. The most 
remarkable feature of our formalism is the greatly reduced uncertainty due to variations of the renormalisation 



scale. We have also implemented our formalism for the case of collisions with hadronic initial states 19 
where we find a similar behaviour. 
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